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INTENSITY MAPPING OF LYMAN-ALPHA EMISSION DURING THE EPOCH OF REIONIZATION 
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ABSTRACT 

We calculate the absolute intensity and anisotropies of the Lyman-a radiation field present during 
the epoch of reionization. We consider emission from both galaxies and the intergalactic medium 
(IGM) and take into account the main contributions to the production of Lyman-a photons: recombi- 
nations, collisions, continuum emission from galaxies and scattering of Lyman-n photons in the IGM. 
We find that the emission from individual galaxies dominates over the IGM with a total Lyman-a 
intensity (times frequency) of about (1.43 — 3.57) x 10 -8 erg s _1 cm~ 2 sr _1 at a redshift of 7. This 
intensity level is low so it is unlikely that the Lyman-a background during reionization can be estab- 
lished by an experiment aiming at an absolute background light measurement. Instead we consider 
Lyman-a intensity mapping with the aim of measuring the anisotropy power spectrum which has rms 
fluctuations at the level of 1 x 10~ 16 [erg s _1 cm~ 2 sr~ x ] 2 at a few Mpc scales. These anisotropies 
could be measured with a spectrometer at near-IR wavelengths from 0.9 to 1.4 /im with fields in the 
order of 0.5 to 1 sq. degrees. We recommend that existing ground-based programs using narrow 
band filters also pursue intensity fluctuations to study statistics on the spatial distribution of faint 
Lyman-a emitters. We also discuss the cross-correlation signal with 21 cm experiments that probe HI 
in the IGM during reionization. A dedicated sub-orbital or space-based Lyman-a intensity mapping 
experiment could provide a viable complimentary approach to probe reionization, when compared to 
21 cm experiments, and is likely within experimental reach. 

Subject headings: cosmology: theory — large scale structure of Universe — diffuse radiation 



1. INTRODUCTION 



The epoch of reionization (EoR) is a crucial stage in 
the history of galaxy formation, signaling the birth of the 
first luminous objects, during which the universe went 
from completely neutr al to almost completely ionized 
(Barkana & Loeb|2001 ). This phase has been largely un- 



explored so far, although current observations suggest it 
was r easonably extended (Komatsu et al. 2011 Fan et al. 
|2006[ ) and a wide variety of observational avenues are be- 
mg explored to probe it. In particular the 21-cm line of 
neutral hydrogen is now understood to be a promising 
tool to study reionization and to understand the forma- 
tion and evolution of galaxies during that epoch (see e.g. 
Furlanetto et al. 2006 ) . It is also now becoming clear that 



we need complimentary data in order to obtain extra in- 
sight into the sources of reionization. Such complimen- 
tary data could also aid in the interpretation of the HI 
signal by allowing ways to pursue cross-correlations and 
providing ways to reduce systematics and foregrounds 
encountered in 21-cm observations. 

Recently, intensity mapping of other atomic and molec- 
ular lines at high redshifts, in particular CO and CII 
( Gong et al.]2012l [20TT| |Lidz et al.|20lH |Visbal fc Loeb 
2U10p , has been proposed as a probe of reionization. In 



anisotropies of the Ly-a line and suggest the latter as 
a new probe of reionization. In particular we suggest 
that it will be possible to measure the amplitude of 
the Ly-a intensity fluctuations with a narrow-band spec- 
trometer either from the ground with a suppression of at- 
mospheric lines or from the orbital/sub-orbital platform. 

The Ly-a line, corresponding to transitions between 
the second and first energy level of the hydrogen atom, 
has a rest wavelength of approximately Xh ya = 1216 A. 
The signal present during reionization is observable in 
near-IR wavelengths today. Existing imaging observa- 
tions made with narrow-band filters on 10m class tele- 
scopes focus on individual galaxy detections and are lim- 
ited to a handful of narrow atmospheric windows at near- 
IR wavelengths. Given the strength o f the line, it ha s 
now been seen in galaxies at z ~ 6.98 ( lye et aL||2006[), 
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this work we study the viability of also using intensity 
mapping of the Lyman-a (Ly-a) line as an additional 
probe. For this study we include several Ly-a emission 
mechanisms involving both individual sources of emission 
such as galaxies and the emission and scattering associ- 
ated with the intergalactic medium (IGM). 
We consider both the integrated intensity and 



narrow-band 
Ly-a emitters have led 
number of galaxies at redshifts 5.7, 6.6, 7.0 and 7.7 to al 
low constraints on the bright-end of the Ly-a lum inosity 
function (L F) and its redshift evolution (e.g. |Ouchi 
Ota et al.||2010[ |Tamguchi et al ~ 
Shibuya et al.||2011|). 
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Observations of the 
F indicate a decrease m the Ly-a intensity from 
redshift 5.7 to 7.0. This would require a strong evolution 
of the Ly-a emitters population, which is not predicted 
by most recent galax y evolution models ( Ota et al.|2010 



Shibuya et al.|[20TT ), or could be explained as the result 



of an increase in the fraction of IGM neutral hydrogen 



2 



that would absor b or scatter Ly-Q photons from th e 
observed galaxies (jHaiman et al.pOOO Ota et aLj [2008 ) . 

The scattering of Ly-a photons by neutral hydrogen m 
the ISM (interstellar medium) and the IGM is expected 
to disperse th e photons in both frequency and direction 
(Santos 2004). Such scattering could considerably de- 
crease the Ly-a intensity per frequency bin from an in- 
dividual galaxy, making the detection of most of the high 
redshift galaxies impossible with current instruments. 
Exact calculations related to scattering are a difficult 
problem to solve analytically and in simulations the scat- 
tering problem requires ray tracing of pho tons through 
the neutral medium in a simulation box ( |Zheng et al 



20101. While scattering makes individual galaxies dim- 



mer, intensity mapping of the Ly-a line at high redshifts 
can be an improvement over the usual experiments that 
make detections of Ly-a emission from point sources and 
are only sensible to the strongest Ly-a emitters. These 
are likely to be some of the brightest star-forming galax- 
ies, however, any dust that is present in such galaxies, 
especially during the late stages of reionization, is likely 
to suppress the Ly-a line. An experiment targeting the 
integrated emission will be able to measure all the sources 
of Ly-a photons in a large region and will be sensitive to 
the extended, low surface brightness Ly-a emission that 
is now kno wn to generally form around star-form ing re- 
gions (e.g., |Steidel et al.pOllj jZheng et al.||2011| . The 
anisotropy power spectrum of Ly-a intensity then would 
be a probe of the Ly-a halos around star- forming galaxies 
present during reionization. The cross-correlation with 
the 21-cm data could provide a direct test on the pres- 
ence of neutral hydrogen in the extended Ly-a halo. 

The paper is organized as follows: in the next section 
wc estimate the contribution to the Ly-a emission from 
galaxies. In section [3] we analyze the contributions to 
the Ly-a emission from the IGM. In section [4] we calcu- 
late the intensity of the Ly-a signal as well as its power 
spectrum using a modified version of the code SimFast21 



(Santos et al 
correlation of 



2010 2011). In section [5] we discuss the 
y-a intensity maps with the 21 cm signal 



and finally in section 6 we comment on the experimental 
feasibility of measuring the Ly-a intensity power spec- 
trum. 

2. LYMAN-a EMISSION FROM GALAXIES 

The observed Ly-a flux is mainly the result of line 
emission from hydrogen recombinations and collisional 
excitations in the interstellar clouds or in the IGM pow- 
ered respectively by UV emission or UV and X-ray emis- 
sion from galaxies. High energy photons emitted by stars 
ionize hydrogen that then recombines to emit a rich spec- 
trum of lines including a Ly-a photon ( Gould & Wein- 



berg 1996 Fernandez fc Komatsu|[2006[ ) . Moreover, the 
electron ejected during this ionization heats the ISM or 
the IGM, increasing the probability o f Ly-a photon emis- 
sion caused by collisional ex citation (Gould & Weinberg 
|1996[[Cantalupo et al.|20 08). There is also a small contri- 
bution to the lyman alpha flux originated in the contin- 
uum emission from stars between the Ly-a line and the 
Lvman- limit ( Chuzhoy fc Zheng||2007| |Barkana fc Loeb 



2005 1 plus Ly-a trom continuum tree- free or tree- bound 



emission as well as 2-photon emission during recombi- 
nations. This continuum will also make contributions 
to a given observation from lower redshifts besides the 



"Ly-a " redshift ( Cooray et al.| 2012 ) which will confuse 
the Ly-a signal. However, due to the smoothness of that 
continuum across frequency, we expect it should be pos- 
sible to remove this contribution, for instance, by fitting 
a smooth polynomial in frequency for each pixel. 

Another source of Ly-a emission in the universe is cool- 
ing of gas that has suffered in-fall into a dark matter halo 
potential well. Several studies show that muc h of this 
cooling is made in the form of Ly-a emission flHaiman] 
eTaL]|2000] I Fardal et al"1[200Tl |Di.jkstra et al.||2ti06a|b[ 
Dayal et al.||2010| |Latif et al ||201ip . Cold gas is used 
by galaxies as fuel to form stars so there is a relation 
between the star formation rate (SFR) of a galaxy and 
the Ly-a flux emitted as gas cools in that galaxy. 

Since emission of Ly-a radiation is closely connected 
with the star formation, the contribution from the sev- 
eral mechanisms by which Ly-a radiation is emitted in 
galaxies and in the IGM can be related to the SFR of in- 
dividual galaxies or galaxy samples. In order to calculate 
the emission of Ly-a radiation from the IGM during the 
EoR we also need to know the ionized fraction of hydro- 
gen as well as the temperature of the gas in the IGM. Un- 
fortunately _j20tli_the^equant^ const r ained 



at z ^ 6 ([Larson et aL 201 1} Ouchi et al 



et al. 201 lp. Since hydrogen ionization shou 



2010| [Zahn 
d be a con- 
sequence of stellar ionization/X-ray emission, we can in 
principle estimate it by following the SFR history and 
making sure that the resulting evolution of hydrogen ion- 
ized fraction is consistent with current constraints on the 
CMB optical depth. 

In order to obtain the SFR of galaxies at the high red- 
shifts during the epoch of reionization we make use of 
parametrizations that reproduce a correct reionization 
history. Our parametrizations are non linear i n a s imilar 
way to the relations found i n the Guo et al. (2011) and 



in the De Lucia fc Blaizot] ( 2007 ) galaxies catalogs de 
rived respectively trom then igh resolution Mille nnium II 



(|Boyla n-Kolc hin et a,l.||20"09| ) and Millennium I ( Springel 
|et al.||2005| simulations. Such relations, when available 
from observations, make an improvement on the models 
instead of relying purely on theoretical calculations and 
semi- numerical simulations to pred i ct all of the observ a- 
tions ( |Mesinger fc Furlanetto"||2007| |Santos et al.|[2010] ). 

There are additional sources of radiation contribut- 
ing to the Ly-a emission, such as a strong non-local 
sources of ionizing photons as expected from quasars, 
which would emit a large amount of energy in X-ray pho- 
tons that would be able to ionize several neutral atoms 
giving origin to a locally strong Ly-a emission from re- 
combinations. However, since the number of quasars is 
very small compared to the number of normal galaxies 
at the redshifts we are considering, we will neglect their 
contribution in the following calculations. 

We encourage future works on Ly-a intensity to see 
if the shape of the power spectrum and other statistics 
can be used to choose between reionization histories that 
involve both galaxies and quasars. 

In the following sub-sections we discuss in more detail 
the four processes for Ly-a emission from galaxies: re- 
combinations, excitations/relaxations, gas cooling, and 
photon emission from continuum processes. 

2.1. Lyman-a emission from hydrogen recombinations 



3 



Assuming ionizing equilibrium, the number of recom- 
binations in galaxies are expected to match the number 
of ionizing photons that are absorbed in the galaxy and 
does not escape into the IGM. Depending on the temper- 
ature and density of the g fraction of the radiation 
due to these recombinations is emitted in the Ly-a line. 

In the interstellar gas, most of the neutral hydrogen 
is in dense clouds with column densities greater than 
3 x 10 18 cm~ 2 . These clouds are optically thick to 
Ly-a radiation and Lyman photons are scattered in the 
galaxy several times before escaping into the IGM. Such 
multiple scatterings increase the probability of absorp- 
tion. Assuming that these clouds are spherical and that 



Gould fc 
re pro ba- 



the gas tempera ture is of the order of 10 4 K 
Weinberg ( 1996 ) used atomic physics to study t 
bility of the Ly-a emission per hydrogen recombination 
They estimated that a fraction f rec ps 66% of the hy- 
drogen recombinations would result in the emission of a 
Ly-a photon and that most of the other recombinations 
would result in two-photon emission. These fractions 
should change with the temperature and the shape of 
the cloud, but such variations are expected to be small. 
Other calculations yield fractions between 62% and 68% 
according to the conditions in the cloud. In this paper we 
have chosen to use a value of f rec = 66% since the overall 
uncertainty on this number is lower than the uncertainty 
on the number of hydrogen recombinations. 

The absorption of Ly-a photons by dust is difficult to 
estimate and ch anges from galaxy to galaxy, | Gould fc] 
Weinberg ( 1996 ) estimated that for a cloud with a col- 
umn density N ~ 10 19 cm -2 , the dust in the galaxy ab- 
sorbs a fraction /dust ~ 4% of the emitted Ly-a photons 
before they reach the galaxy virial radius however recent 
observations of high redshift galaxies indicate a much 
higher /dust- In this study we will use a redshift param- 
eterization for the fraction of Ly-a photons that are not 
absorbed by dust fhya — l~/dus t that is doub l e the value 



predicted by the study made by |Hayes et al 

fhya(z) 

where C dust 



(2011) 



Cdust x l(T d (l 
3.34 and £ = 2.57. 



z) 



(1) 



The |Hayes et al. 
(2011) parameterization was made so that fhya gives 



the difference between observed Ly-a luminosities and 
Ly-a luminosities scaled from star formation rates as- 
suming that the Ly-a alpha photons emitted in galaxies 
are only originated in recombinations. The high redshift 
observations used to estimate fhya are only of massive 
stars while the bulk of Ly-a emission is originated in the 
low mass stars that cannot be det ected by current sur- 
veys. According to several studies ( Forero- Romero et al. 
20111, fhya decreases with halo m ass, so it is possible 



that 



(2011) 



it is being underestimated in Hayes et al 
which is why we decided to use a higher fhya ■ Uur results 
can however be easily scaled to other /dust evolutions. 
The number of Ly-a photons emitted in a galaxy per 



second, Ni 



Lyaj 



that reach its virial radius is therefore 



given by 

Nhya = Affe/rec X fhya X (1 - / osc ) X iV ion , (2) 

where A^c — accounts for the fraction of photons 

that go into the ionization of helium (Y p is the mass 
fraction of helium), 7V; on is the rate of ionizing photons 



z a 




P 


/ eBC (M = 10 10 AT©) 


10.4 2.78 x 10" 


-1 


0.105 


0.732 


8.2 1.30 x 10" 


-2 


0.179 


0.449 


6.7 5.18 x 10- 


-3 


0.244 


0.240 


5.7 3.42 x 10" 


-3 


0.262 


0.240 



TABLE 1 

Fits to the escape fraction of UV radiation from galaxies 
as a function of redshift (based on razoumov & 
Sommer-Larsen 2010). 

emitted by the stars in the galaxy and / osc is the fraction 
of ionizing photons that escape the galaxy into the IGM. 

The ionizing photon escape fraction depends on con- 
ditions inside each galaxy and is difficult to estimate, 
especially at high redshifts. The precise determination 
of its value is one of the major goals of future observa- 
tions of high redshift galaxies at z > 7. This parameter 
can be measured from deep imaging observations or can 
be estimated from the equivalent widths of the hydrogen 
and helium balmer lines. The ionizing photon escape 
fraction dependence with the galaxy mass and the star 
formation rate, as a function of redshift, has been esti- 
mated using simulations that make several assumptions 
about the intensity of this radiation and its absorption 
in the interstellar medium. However, for the halo virial 
mass range, 10 8 M Q to 1O 13 M , and during the broad 
redshift range related to the epoch of reionization, there 
are no simulations that cover the full parameter space. 
Moreover the limited simulations that exist do not al- 



ways agree with each other dGnedin et , al.|2008 Wise & 



CeHl[2009l |Fernande z-Soto et al. 200 31 ISiana et al ||2007t 
Haardt fc Madaul|2012p .~Razoumov fc Sommer-Larsen 



] 2010 1 computed the escape traction ot U V radiation tor 
the redshift interval z — 4 to z = 10 and for halos of 
masses from 10 78 to 10 115 M Q using a high-resolution 
set of galaxies. Their simulations cover most of the pa- 
rameter space needed for reionization related calculations 
and their escape fraction parameterization is compatible 
with most of the current observational results. Thus, we 
use it for our calculations here. 



According to Razoumov & Sommer-Larsen (2010) sim- 
ulations, the escape traction of ionizing radiation can be 
parameterized as: 



/ e8c (M,z)=exp \-a(z)M^ 



(3) 



where M is the halo mass, a and (3 are functions of red- 
shift (Table [l]). The number of ionizing photons emit- 
ted by the stars in a galaxy depends on its star forma- 
tion rate, metallicity and the stellar initial mass function 
(IMF). Making reasonable assumptions for these quanti- 
ties we will now estimate iV; on . Since this UV emission is 
dominated by massive, short lived stars, we can assume 
that the intensity of ionizing photons emitted by a galaxy 
is proportional to its star formation rate. In terms of the 
star formation rate in one galaxy, 



Ni on — Qi on X SFR, 



(4) 



where Qi on is the average number of ionizing photons 
emitted per solar mass of star formation. This can be 
calculated through: 



/^ max *(M)Q+(M)U{M)dM 



Qic 



*(M)MdM 



(5) 



4 



where *(M) = KM~ a is the stellar IMF, if is a constant 
normalization factor and a is the slope of the IMF. In 
our calculation we used a Salpeter IMF, with a = 2.35. 
t±(M) is the star lifetime and Q+(M) its number of ion- 
izing photons emitted per unit time. The values of 
an d t+ were calcul ated with the ionizing fluxes obtained 
by Schaerer ( 2002 ) using realistic models of stellar popu- 
lations and non-LlE atmospheric models, appropriated 
for POP II stars with aZ t = 0.02Z© metallicity. 

Assuming that ionizing photons are only emitted by 
massive OB stars sets a low mass effective limit for the 
mass of stars contributing to the UV radiation held of a 
galaxy. This limit is a necessary condition for the star 
to be able to produce a significant number of ionizing 
photons. For the stell ar population used for this wor k 
we take M min w 7 M ( lScnaerer|2002]|Shurl et al.|2011j ). 
The integration upper limit is taken to be iW max = 150 
M Q . In this paper we calcula ted Q,; nT] , using t he parame- 
terization values published in Schaerer (2002 ). The num- 
ber of ionizing photons per second emitted by a star as 
a function of its mass is given by: 

logi [Q*/s _1 ] = 27.80 +30.68:c - 14.80x 2 (6) 
+ 2.5x 3 for 7M < M* < 150 M 

where x = \og w (M ir / M & ) and the star's lifetime in years 
is given by: 



log 10 [**/yr] = 9.59 - 2.79.x + 0.63a; 2 



(7) 



The use of these parameter s results in Q lon « 5.38 x 10 60 



m; 



In 



Shull et al. 



of a different model 
(R. S 



(2011 ) it has been suggested the use 
for stellar atmosphere and evolution 
Sutherland & J. M. Shull, unpublished) which 



yields Q ion 3.97 x 10 60 M~ . This may imply that 
the stellar emissivity we calculated is an overestimation 
and that consequently our Ly-a flux powered by stellar 
emission may be overestimated by about 35%. This is 
comparable to other large uncertainties, such as the ones 
in the parameters / esc and /dust ■ The Ly-a luminosity is 
calculated assuming that the Ly-a photons are emitted 
at the Ly-a rest frequency, vq = 2.47 x 10 15 Hz with an 
energy of E\, ya = 1.637 x 10~ n erg. To proceed, we will 
assume that the SFR for a given galaxy is only a func- 
tion of redshift and the mass of the dark halo associated 
with that galaxy. The Ly-a luminosity due to recombi- 
nations in the interstellar medium, L^ L , can then be 
parameterized as a function of halo mass and redshift as 



^rec (^"i z ) — EhyaNhya 



(8) 

: 1.55 x 10 42 [1 - f esc (M, z)] / Lya (z) S A F . R(M 'f erg s~K 

M yr- 1 



2.2. Lyman-a emission from excitations/relaxations 

The kinetic energy of the electron ejected during the 
hydrogen ionization heats the gas and assuming ther- 
mal equilibriu m this heat is emitted as r adiation. Using 
atomic physics, Gould & Weinberg (1996) estimated that 
for a cloud with an hydrogen column density of rs 10 
cm~ 2 , the energy emitted in the form of Ly-a photons 
is about 60% for ionizing photons with energy E v 
E v < AE V 



< 

and w 50% for photons with energy E v > 
= 13.6 eV is the Rydberg energy. 
The remaining of the energy is emitted in other lines. 



4:E Vlim , where E Vl 



Using the spectral energy distribution (SED) of galax- 
ies with a meta llicity Z — 0.02Z Q from the code of 
Maraston (2005) we estimated that the average ioniz- 
ing photon energy is E v = 21.4 eV and that more than 
99% of the phot ons have an energy lower than AE vliTrL . 
According to the Gould & Weinberg ( 1996 ) calculation, 
the fraction of energy of the U V photon that is emitted 
as Lyman alpha radiation due to the collisional excita- 
tions/relaxations is given by: 



E c * c /E v ~ 0.08 + 0.1 1 



2^1ir 



0.1 



(9) 



For a cloud with the properties considered here this yields 
an energy in Ly-a per ionizing photon of i? cxc 2.14 eV 
or 3.43 x 10 -12 erg. This results in an average of 0.16 
Ly-a photons per ionizing photon. 
Finally, the Ly-a luminosity due to excitations in the 



ISM, L^ce , is then 
-GAL 



L^ L (M,z) = [l-f eBC (M,z)]f hya (z)A 



He 



(10) 



4.03 xlO 41 [1- f esc (M,z)]f Lya (z) 



SFR(M, z 



M yr 



j- erg s 



where again it is assumed to be a function of the star 
formation rate. 

2.3. Lyman-a emission from gas cooling 

During the formation of galaxies, gas from the IGM 
falls into potential wells composed mainly by dark matter 
which collapsed under its own gravity. The increase in 
the gas density leads to a high rate of atomic collisions 
that heats the gas to a high temp erature. According 
to the study of 
potential wells t 



Fardal et al. (2001) most of the gas in 
lat collapses under its own gravity never 
reaches its virial temperature and so a large fraction of 
the potential energy is released by line emission induced 
by collisions and excitations from gas with temperatures 
Tjc < 2 x 10 4 K. At this temperature approximately 50% 
of the e nergy is emit t ed in Ly-a alone. 

From Fardal et al. (2001) we can relate the luminosity 
at the Ly-a frequency due to the cooling in galaxies to 
their baryonic cold mass, 



tcool' usm S : 



log 10 (^ L ) 



1.521og 10 (M, 



bar ^ 
cool I 



26.32, (11) 



where both the luminosity and the mass are in solar 
units. To relate this baryonic cold mass to a quantity we 
can use in our models, we used the relation between cold 
baryonic mass an d the halo mass from the galaxies in the 
Guo et al. (2011) catalog. From the equation above, we 
can then obtain an expression for the luminosity, which 
can be fitted by: 



L^i(M) w 1.69 x 10 ao h ya (z) [ 1 + 



35 . 



M 

w 



(12) 



x 1 + 



M 
2 x 10 10 



2.1 



1 



M 
3 x 10 1 



with M in units of M Q . The relation between the cold 
gas mass and the mass of the halo shows very little evo- 
lution with redshift during reionization. Thus we expect 
the relation in equation 13 to only depend on redshift 



due to the redshift evolution of fhya- 



5 



2.4. Contributions from continuum emission 

Continuum emission can also contribute to the 
Ly-a observations. These include stellar emission, free- 
free emission, free-bound emission and two photon emis- 
sion. Photons emitted with frequencies close to the 
Ly-n lines should scatter within the ISM and eventually 
get re-emitted out of the galaxy as Ly-a photons. Other- 
wise they will escape the ISM before redshifting into one 
of the Ly-n lines and being reabsorbed by a hydrogen 
atom. 

The fraction of photons that scatter in the galaxy can 
be estimated from the intrinsic width of t he Lyman alpha 
line which has w 4A ( Jensen et al. 1|2012[ ). We calculated 



the stellar contribution assuming an emission spectrum 
for stars with a metallicity of Z+ = 0.02Z Q estimated 
with the code from Maraston ( 2005 ) that can be approx- 
imated by the emission of a black body with a temper- 
ature of 6.0 x 10 4 K for hv < 13.6eV. The number of 
stellar origin Ly-a photons per solar mass in star forma- 
tion obtained with this method is: 



^^ = 4.307 



dv 



, " Mo 1 (13) 



= 9.92 x 10 58 Mq 1 . 



We note that we are not accounting for the higher opacity 
at the center of the Ly-a line which should push the pho- 
tons out of the line center before exiting the star and so 
we may be overestimating the stellar Ly-a photon emis- 
sion. 

Free-bound emission and free-free emission are respec- 
tively originated when free electrons scatter off ions with 
or without being captured. F ollowing the approach of 
(Fernandez & Komatsu 2006), the free- free and free- 
bound continuum luminosity can be obtained using: 



L v (M,z) = V sphere (M,z)e u 



(14) 



where V sp here is the volume of the Stromgren sphere 
which can be roughly estimated using the ratio between 
the number of ionizing photons emitted and the number 
density of recombinations in the ionized volume, 



V sphere (M,z) = 



Qi on SFR(M, z)(l-/ e8C ) 
n e n p ap 



(15) 



e v is the total volume emissivity of free-free and free- 
bound emission, n p is the number density of protons (ion- 
ized atoms) an d Uj is the case A or case B recombination 
coefficient (see Furlanetto et al. (2006)) 



The volu me emissivity estimated by ( Dopita & Suther- 
Iandl|2003 ) is given by: 



47rn e n P 7 c - 



-hv/kT K 



T. 



1/2 



-Jem 3 s 1 Hz \ 



(16) 



A' 



where j c is the continuum emission coefficient including 
free-free and free-bound emission given in SI units by: 



7c 



5.44 x 10~ 46 



9ff + ^=nl 



-9fb{n) 



(17) 



In here, x n = Ry / {kBTxn 2 ) {k,B is the Boltzmann con- 
stant, n is the level to which the electron recombines 
to and Ry — 13.6eV is the Rydberg unit of energy), 



#// ~ 1.1 — 1.2 and gfb(n) ~ 1.05—1.09 are the thermally 
avera ge Gaunt factors for fre e-free and free-bound emis- 



sion (Karzas & Latter 1961 values from). The initial 



level n is determined by the emitted photon frequency 
and satisfies the condition cR^jn^ < v < cR^/ln — l) 2 
where R^ = 1.1 x 10 7 m _1 is the Rydberg constant. 

The continuum luminosity per frequency inter- 
val {L u ) is related to the Ly-a luminosity emitted 

from the galaxies by: L cont — L u x dv(4A) — 
fL V aQLyaE L y a SFR(M, z), where Qhya is the number 
of emitted lyman alpha photons per solar mass in star 
formation. We then obtain Q^~ ilee = 2.13 x 10 53 Mq 1 

for free-free emission and QLyr b ° Und = 2 - 22 x 1055 M o' 
for free-bound emission. 

During recombination there is also the probability of 
two photon emission and although this photons have fre- 
quencies below the lyman alpha frequency there is a small 
fraction of them of Q^^ hot ° n ^ na ^ are em itted so close 
to the lyman alpha line that are included in the lyman 
alpha intrinsic width. 

The number of lyman alpha photons that can be orig- 
inated due to two photon emission during recombination 
is given by: 



-j2 — photon 
-Lya 



l^Lya 



P{v/v Lya )dv, 



(18) 



where P(y)dy is the normalized probability that in a two 
photon decay one of them is the range dy — dv jvi ya and 
1 — fiya ~ 1/3 is the probability of 2 photon emission 
during an hydrogen n=2— >1 tran sition. The probability 
of two photon decay was fitted by |Fern andcz fc Komatsu| 
( |2006| using Table 4 of Brown & Mathews (1970) as: 

P(y) = 1.307 - 2.627(y - 0.5) 2 + 2.563(y - 0.5) 4 (19) 
-51.69(y- 0.5) 6 

Finally, the different contributions to the total 
Ly-a luminosity from galaxies due to continuum emis- 

ainn T GAL _ r stellar i r fi 
sion, l^ cont — i^ cont 

are given by: 



- free— free, r free— bound , r 2— photon 
^cont "'"-^cont "r-^cont 5 



L±t ar (M, Z) = / L yaO L ya ar ^LyaSFR(M, z) 



(20) 



M© yr- 



for stellar emission, 



C" « '-) = /L»„<?f y T'EL,„SFR(M, z) (21) 



M Q yr- 



for free-free emission, 



- free— bound / 



(M, Z) = /LyaQLyr£LyaSFR(M, z) (22) 



:1.47 X 10 37 /Lya 



SFR(M, z) , 
M yr 1 



for free-bound emission and 



- 2 — photon 



(M, Z) = / L y a g L t y °a ar ^LyaSFR(M, z) 



(23) 



, oo , SFR(M, z) , 
' 2.41 x 10 38 /Lya M 1 I/ erg s' 1 
M yr 1 



G 



for 2-photon emission. 

Note that here we are only considering the part of the 
continuum emission from galaxies that could contribute 
to the same "Ly-a redshift". There will be a continuum 
emission spectrum with frequencies below the Ly-a line 
from the mechanisms above that will contribute to the 
same observation from lower redshifts and will generate 
a "foreground" to the Ly-a signal that needs to be re- 
moved. This should be possible due to the smoothness 
of this background across frequency in the same man- 
ner as foreground s of the 21-cm signal are removed (e.g. 
Wang et al.|[2006| ). 



2.5. Modeling the relation between star formation rate 
and halo mass 

Simulations of galaxy formation and observations indi- 
cate that the star formation of a halo increases strongly 
for small halo masses but at high h alo masses (M > 
lO^Mp)) it becomes almost c onstant ( Conroy & Wech- 
s1ct|[2009| [Popesso et al.|[2012] ). 

In order to better estimate and constrain the SFR of 
a halo we used three non linear SFR versus Halo Mass 
parameterizations that are in good agreement with dif- 
ferent observational constraints. In Siml we adjusted 
the SFR to reproduce a reasonable reionization history 
and a Ly-a Luminosity Function evolution compatible 
with different observational constraints, in Sim2 we ad- 
justed the SFR vs halo mass rela tion to the parameteri- 
zations from the Guo ct al. ( 201 1 j) galaxies catalo gue (low 



halo masses) and the De Lucia & Ulaizot p007 ) galaxies 
catalogue (high halo masses). t>im2 results in an early 
reionization history with an optical depth to reionization 
compatible with the low bound of the current observa- 
tional constraints. Finally SimS has the same halo mass 
dependence as Sim2 but evolves with re dshift in a simi- 
lar way to th e De Lucia & Blaizot ( 2007 ) and to the Guo 
et al. (2011) galaxy catalogues. 

We parametrized the relations between the SFR and 
halo mass as: 

SFR ^. (2 .8xio-)ivrx 

M /yr V ' 



where a = 2.8, b = — u.y4, a = 
and c 2 = 7 x 10 10 M Q for Siml 



(24) 

= -0.94, d = -1.7 , ci = 1 x 10 9 M Q 

' M„ fnr Rim 1 



M Q /yr 

' M\ 6 / 



1 



ci 



C2 



AT 

c 3 



(25) 

where a = 2.59, b = -0.62, d = 0.4, e = -2.25, a = 
8 x 10 8 M , c 2 = 7x 10 9 Mq and c 3 = 1 x 10 11 M for 
Sim2 and 



SFR(M,,) =225 xl0 _ 2 Q75x 
M /yr 



(z - 7)) M a x 
L + f)', (26) 



where a = 2.59, b = -0.62, d = 0.4, 
8 x 10 8 M Q , c 2 = 7 x 10 9 M and c 3 = 
Sim3. 

Figure [l] shows these relations. 
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Fig. 1. — Star formation rate ver sus halo mass. Th e dotted lines 
show the relations taken from the |Guo et al.| ( |201 1 [ > catalogue for 
low halo masses at z = 6 (bottom dotted linej and z = 8 (up- 
per dotted line), the yellow crosses show the relation taken from 
the DeLucia catalogue for high halo masses at z = 7. The dash- 
dotted, solid and dashed lines show the parameterizations used in 
simulations Siml, Sim2 and SimS respectively for z = 7. 




Fig. 2. — Star formation rate density evolution as a function of 
redshift. The blue solid line, the green dashed line and the black 
dashed dotted line were obtained from simulations made using the 
SimFast21 code (for informations about the code see section |4| and 
Santos ct al. 2010) and the SFR vs. halo mass relations rrom 
cquations |24| [25] and 1 2 6 1 The red dots are observational constraints 
derived from the UV luminosities corrected for dust extinction from 
|Bouwens et a l. (2011c). Please note that these observational values 
correspond to hign mass galaxies while our results integrate over 
the halo mass function starting at ~ 10 8 solar masses (which at 
redshift 7 corresponds to star formation rates of 6.41 X 10 — 5 Mq 
s- 1 for Siml, 7.83 xl0~ 6 Mq s" 1 for Sim2 and 1.1 xlO" 5 Mq 
s — 1 for SimS), so our star formation rate densities are expected to 
be higher. 
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In figure [2j the strong decline in the observational 
SFRD from z m 8 to z w 10, imposed by the obser- 
vational point at z — 10.3, was obtained with the ob- 
servation of a single galaxy using the Hubble Deep Field 
2009 two years data ( jBouwens et al.||20lia| |Oesch et al. 
2012 1. It was argued m|Bouwcns ^!TaE|(|2011b[ ), based on 
an analytical calculation, that even with, such low SFRD 
at high redshifts it was possible to obtain an optical 
depth to reionization compatible with the value obtained 
by WMAP (t = 0.088 ± 0.015) flKomatsu et"aJT|[25TH >. 
However, this derivation would imply a high escape frac- 
tion of ionizing radiation and that reionization would 
end at z ~ 8 which is hard to reconcile with the con- 



straints from observations of quasars spe ctra (Mesinger 
fc Haiman||2007| |Zaldarriaga et al.||2008| . Our SFHDs 
are considerably higher than the current observational 
constraints, although the difference can be explained by 
a systematic underestimation of the SFR in observed 
galaxies. Moreover, current observations only probe the 
high mass end of the high redshift galaxies mass func- 
tion which will underestimate the SFR density (also the 
obtained SFRs have very high error bars due to uncer- 
tainties in the correction due to dust extinction, the red- 
shift and the galaxy type). In the following sections the 
results shown were obtained using Siml unless stated 
otherwise. 

2.6. Total Lyman-a luminosity: comparison with 
observations 

In the previous sections we calculated the 
Ly-a luminosity as a function of the SFR for sev- 
eral effects. The commonly us ed "empirical" re lation 
between these two quantities is ( Jiang et al.||20li| ) 



L Gal = 1.1 x 10 4 



, SFR(M, z) 
Moyr- 1 



erg s 



(27) 



and it is based on the relation betwee n SFR and the 
Ha luminosity from Kennicutt ( 1998a[) and in the line 
emission ratio of Ly-a to Ha in case B recombinations 
calculated assuming a gas temperature of 10 4 K. This 
empirical relation gives the Ly-a luminosity without dust 
absorption (we have labeled it K98 for the remainder of 
the paper). 

Our relation between luminosity and star formation is 
mass dependent (both from the escape fraction as well as 
due to the expression from the cooling mechanism), so in 
order to compare it with the result above, we calculate: 



A(z) 



(L Gal (M,z)) 
(SFR(M, z)) ' 



(28) 



where the average (x) of quantity x is done over the halo 
mass function for the mass range considered. The results 
are presented in table [2] for a few redshifts. 

Although our Ly-a luminosities per SFR are slightly 
higher, at least for low redshifts, we point out that the 
"empirical" relation is based on a theoretical calculation 
that only accounts for Ly-a emission due to recombina- 
tions. Moreover the observational measurements of Ha 
and Ly-a are primarily made at low redshifts, where the 
absorption of Ly-a photons by dust in galaxies is ex- 
pected to be high. Our relation has the advantage of 
evolving with redshift since it accounts for the evolution 
of the escape fraction of ionizing photons and for the 



A rcc (z) A CX c(z) 

4.4 x 10 41 1.1 x 10 41 



Acool(z) A con t(z) Atotal(z) 



1.2 x 10 42 3.2 x 10 

9.3 x 10 41 
8.5 x 10 41 



,41 

2.4 x 10 41 
2.2 x 10 41 



1.3 x 10 39 

7.8 x 10 38 

4.9 x 10 38 
3.3 x 10 38 



8.1 x 10 40 6.4 x 10 41 

6.4 x 10 40 1.6 x 10 42 

4.9 x 10 40 1.2 x 10 42 

3.6 x 10 40 1.1 x 10 42 



TABLE 2 

Average luminosity per star formation rate (in units 

erg S -1 /M0 yr _1 ) AVERAGED OVER THE HALO MASS FUNCTION 
FOR REDSHIFTS 10. 9, 8 AND 7, FROM TOP TO BOTTOM. 



evolution of the escape fraction of Ly-a photons. This 
z-dependence is not present in the standard empirical re- 
lation. This redshift evolution of the UV photons escape 
fraction is a consequence of the increase in the number 
of massive galaxies with more clumpy structure as the 
redshift decreases. The star forming regions of massive 
galaxies are embedded in clumps and therefore it be- 
comes more difficult for t he ionizing photons to escape 



from such dense regions ( Razoumov & Sommer-Larsen 
20101 |Yajima et al |[20lT| ). The redshift evolution of the 
relation presented in equation 28]justines why a theoret- 
ical calibration between Ly-a luminosity and the SFR of 
a galaxy is useful for our work. 

To check the consistency between our theoretical esti- 
mation of the Ly-a luminosity and the existing observa- 
tions during reionization we show in figure [3] the lumi- 
nosity function (LF) using two of the star formation rate 
vs. halo mass parameterizations presented in section [275] 
This prediction is then compared to Ly-a luminosity 
functions of photo metric identified objects in[Shimasaku 
et al. ( 2006 ) and in Kashikawa et al. ( 2006 ) near the end 
of the reionization epoch. 




log 10 [L Lya >(erg s')] 

Fig. 3. — Ly-a luminosity functions obtained with our calcula- 
tions are shown for redshifts z = 6 (dashed lines) and 2 = 7 (solid 
lines) for Siml (black thick lines) and Sim2 (blue thin lines). The 
green and red circles show the intrinsic (i.e. not affe cted by the 



IGM) Ly-a LF from photometric iden tified objects in |Shimasaku| 
|et al.| \2006\ and in |Kashikawa et al.| | |2006[ l for z = 5.7 and 6.6 
respectively. 



Our luminosity functions were calculated assuming a 
minimum halo mass of 8 x lO 8 Af0 which corresponds to 



a minimum luminosity of 3.72 x 10 erg s for Siml, 



4.49 x 10 ib erg s _i for Sim2 and 6.22 x 10 ib erg s" 1 for 
Sim3. The agreement between our LFs and observations 
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is reasonable for Siml however our Sim2 overpredicts 
the abundance of high luminosity Ly-a emitters. This 
difference can be due to sample variance or a result of 
the high sensitivity of theoretical predictions to several 
parameters in our model. We point out that the lumi- 
nosity range relevant for this comparison falls in a halo 
mass range outside the one for which the escape fraction 
of UV radiation we are using was estimated, so we could 
easily get a better fit between observations and Sim2 
by reducing this escape fraction for high halo masses. 
This difference could also be related with the choice of 
halo mass function. Here we choose the Sheth-Tormen 



halo mass function (Sheth & Tormen 1999) which has 
been shown to fit low-redshitt simulations more accu- 
rately, but it is yet to be established the extent to which 
such a halo mass function can reproduce the halo distri- 
bution during reionization. Other possible explanation 
for this difference is the existence of a small amount of 
neutral gas in the IGM which would severely decrease 
the observed Ly-a luminosity from galaxies. Also, we 
could have decreased the high luminosity end of our lu- 
minosity functions if we had use an Ly-a escape frac- 
tio n that decreased with halo m ass such as the one used 
in ( |Forero-Romero et al. 2011). We do not consider a 
model fit to the data to optimize various parameters in 
our model given that the current constraints on the ob- 
served Ly-a LFs have large overall uncertainties, espe- 
cially considering variations from one survey to another. 

2.7. Lyman-a Average Intensity 

In this section and the next one we will attempt to esti- 
mate the intensity and power spectrum of the Ly-a signal 
using an analytical model. In section [4] we will improve 
the estimation by doing the same calculation using a 
semi-numerical simulation. 

The total intensity of Ly-a emission can be obtained 
from the combined luminosity of Ly-a photons associ- 
ated with different mechanisms described in the previous 
sub-sections, such that: 



'Gal 



CO 



dM d^ L G ^M,z) y[z)Dl 



M„ 



dM 



4nDl 



(29) 



where dn/ dM is the halo mass function (Sheth & Tor 

men 11 QQQh A/f ;o +Vio ViqI/^ maoc A/f — I I) 1:1 

M 



is the halo mass, M max = IO^Mq, 
Djj is the proper luminosity distance 

Fi- 



1999b, M 

mm = -^OB , 

and Da the comoving angular diameter distance 



nally, y(z)_ = d\/dv = AL yQ (l + z) 2 /H(z), where x is 
the comoving distance, v is the observed frequency and 
ALya = 2.46 x 10~ 15 m is the rest-frame wavelength of 
the Ly-a line. 

The evolution of the Ly-a intensity predicted by this 
calculation is shown in figure [3] together with the scal- 
ing expected under the "empirical" relation from [Ken^ 



nicutt 



(1998a) combined with an assumption related to 
the gas temperature. The intensities of Ly-a emission 
from different sources are presented in table [3] for several 
redshifts. 

These intensities can be extrapolated to other SFRDs, 
assuming that the only change is in the amplitude of the 
SFR halo mass relations presented in figure [T] by using 
the coefficients in table |4j 

The intensities from emission at z ~ 7, 8 and 10 are 
9.5 x 10~ 9 , 3.5 x 10~ 9 and 3.2 x 10~ 10 erg s" 1 cm" 2 sr" 1 , 




Fig. 4. — Ly-a Intensity from galaxies in erg s cm sr as 
a function of rcdshift. The black dashed dotted line and the blue 
solid line were obtained using our theoretical calculation of the 
Ly-a luminosity and the SFR halo mass relation from Siml and 
Sim2 respectively. The orange dotted line uses the Ly-a luminosity 
SFR relation based on th e relation between SFR and the Ha lu- 
minosity from |Kennicutt| ( |1998a[ | and the line emission ratio of 
Ly-a to Ha in case B recombinations calculated assuming a gas 
temperature of 10000K (labeled as the K98 relation). The K98 
line is not corrected for dust absorption. 



Source of emission in 

r -1 - 2 -li 
[erg s cm sr J 


i/I„(z=7) 




) 


i/I„(z=10) 


Recombinations 


7.3 x 10" 9 


2.5 X 10" 


-9 


2.3 x 10" 10 


Excitations 


1.9 x icr 9 


6.5 x 10- 


10 


6.0 x lO" 11 


Cooling 


2.8 x icr 12 


1.5 x 10" 


12 


4.7 x 10~ 13 


Continuum 


3.1 x 10~ 10 


3.5 x 10- 


10 


3.0 x 10" 11 


Total 


9.5 x 10" 9 


3.5 x 10" 


-9 


3.2 x 10- 10 



TABLE 3 

Surface brightness (in observed frequency times intensity) 
of Ly-o emission from the different sources in galaxies at 



ar ss 7, 


Z ~ 8 AND Z £3 


10 FOR Siml. 


Redshift 


AsFRD(Siml) 


A SFRD (Sim2) 


10 


1.05 x 10" 03 


2.20 x 10" 03 


9 


2.17 x lO" 03 


3.75 x 10-° 3 


8 


2.34 x 10" 03 


3.15 x 10" 03 


7 


3.42 x 10" 03 


3.42 x 10" 03 



TABLE 4 

Average LY-a Intensity from galaxies per SFRD (As FIiD ) 

IN UNITS erg s- 1 cm- 2 sr~ 1 /AfQj/r -1 , CALCULATED USING THE 
STAR FORMATION RATE HALO MASS RELATION FROM SIMULATIONS 
Siml AND Sim2. 



respectively. Such an intensity is substantially smaller 
than the background intensity of integrated emission 
from all galaxies (around 1 x 10 -5 erg s _1 cm _2 sr _1 



(Madau & Pozzetti 2000|, or from the total emission 
of galaxies during reionization, estimated to be at most 
1 x lO" 6 erg s^cm" 2 sr" 1 (Yan et al. 2012). 

2.8. Lyman-a intensity power spectrum 

The Ly-a emission from galaxies will naturally trace 
the underlying cosmic matter density field so we can 
write the Ly-a line intensity fluctuations due to galaxy 
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Fig. 5. — Bias between dark matter fluctuations and Ly-a surface 
brightness (vl) from galaxies as a function of the galaxy Lyman 
alpha luminosity at redshifts 7, 8 ,9 and 10. 




clustering as 



WgAL = 6 L ya/GAL<5(x), 



(30) 



where /gal is the mean intensity of the Ly-a emission 
line, 5(x) is the matter over-density at the location x, 
and &Lya is the average galaxy bias weig hted by the 
Ly-a luminosity (see e.g. Gong et al.||2011 ). 
Using one of t he re lations between the SFR and halo 
we can calculate th e luminosity and 
pha bias following Visbal & Loeb 



mass from section [275 
obtain the Lyman a 
pOlOb: 



M„ 

A-f m 



dM^L GAL b(z,M) 



dM^L GAL 



(31) 



where b(z. M) is the halo bias and dnj dM the halo mass 
function QSheth & Tormen| 119991 ). We take M min = 
10 8 M Q /ft, and M max = W^M Q /h. The bias between 
dark matter fluctuation and the Ly-a luminosity, as can 
be seen in figure [5j is dominated by the galaxies with low 
Ly-a luminosity independently of the redshift. 

We can then obtain the clustering power spectrum of 
Ly-a emission as 



(32) 



where Pss(z, k) is the matter power spectrum. The shot 
noise power spectrum, due to discretization of the galax 
ies, is also considered here. It can be written as (Gong 
et^[20lT ) 



P>shot 
I 



Ly. 



f Mlnin d M^-(^y( Z )D 
7 Mmax dM\4nDl yK ' 



(33) 



The resulting power spectrum of Ly-a emission from 
galaxies is presented in figure [6| At all scales presented 
the Ly-a intensity and fluctuations are dominated by the 
recombination emission from galaxies. 

3. LYMAN-a EMISSION FROM THE IGM 

The Ly-a emission from the IGM is mostly originated 
in recombinations and collisions powered by the ionizing 
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Fig. 6. — Clustering power spectrum of the Ly-o surface bright- 
ness (vl) from Galaxies at redshifts 7 to 10 (from top to bottom), 
from several sources: collisions and excitations, recombinations and 
continuum emission with frequencies inside the Ly-a width. The 
power spectra from cooling emission is not shown since it is sev- 
eral orders of magnitude smaller than the contributions from the 
other sources. Also shown are the total power spectra (clustering 
(solid black line) and shot noise power spectra (dotted black line)) 
of the total contribution for Ly-a emission in galaxies predicted by 
our theoretical calculation and total Ly-ct clustering power spectra 
predicted using the ftf98 relation (orange solid line). 



background. These processes are similar to the ones de- 
scribed inside the galaxies, although, since the physical 
conditions of the gas in the IGM are different from those 
in the ISM, the intensity of Ly-a emission can no longer 
be connected to the ionizing photon intensity using the 
previous relations. The biggest challenge in doing these 
calculations is to connect the IGM ionizations and heat- 
ing of the gas to the emission of ionizing radiation and 
the star formation rate assumed in the previous sections. 
Moreover, in the IGM, we also have to take into account 
the contribution of continuum radiation from stars be- 
tween the Ly-a and the Lyman limit which redshifts into 
the Ly-a line. 

In a schematic view, we have to take into account the 
following processes, 

1. The amount of energy in UV photons that escapes 
the galaxy 
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2. This energy will then be distributed in the IGM 
into: 

(a) ionizations 

(b) direct excitations (followed by emission, par- 
tially into the Ly-a line) 

(c) heating of the gas 

3. Taking into account the state of the IGM in terms 
of temperature and ionization, we can then further 
determine how much it will radiate through the 
Ly-a line from: 

(a) Recombinations 

(b) Radiative cooling (usually through excitations 
followed by decay in several lines including 
Ly-a ) 

4. The amount of Ly-n photons that escape the 
galaxy, re-scattering in the IGM into Ly-a photons 

The proper calculation of all these processes will require 
simulations which we will address in section |4j In the fol- 
lowing sub-sections we review the contributions through 
analytical calculations in order to get a better under- 
standing of the dominating effects. 

3.1. Lyman-a emission from hydrogen recombinations 

The UV radiation that escapes the interstellar medium 
into the intergalactic medium ionizes low density clouds 
of neutral gas. Part of the gas in these clouds then re- 
combines giving rise to Ly-a emission. The radiation 
emitted in th e IGM is often referred to as fluorescence 
( |Santos|[2004| . The comoving number density of recom- 
binations per second in a given region, n rec , is given by: 

n rec (z) = an e (z)n HU (z) (34) 

where a changes between the case A and the case B re- 
combinations coefficient, tihii = x iT~jjj^ is the ionized 
hydrogen comoving number density (xi is the ionization 
fraction, n& the baryonic comoving number density). The 
free electron density can be approximated by n e — Xin^. 

The recombination coefficients are a function of the 
IGM temperature, Tx- The case A comoving recombi- 
nation coefficient is appropriate for the hig hly ionized 



the pa rameterization for f rec made by [Cantalupo et al 



low redshift Universe |Furlanetto et al. ( 2006[ ), 



a A « 4.2 x 10- 13 (r K /10 4 K)-°- 7 (l + zf 



cm 3 s 



(35) 



while the case B comoving recombination coefficient is 
appropriate for the high redshift Universe. 

a B « 2.6 x 10" 13 (T K /10 4 K)"°- 7 (1 + zf cm 3 s _1 . (36) 

The use of a larger recombination coefficient when the 
process of hydrogen recombination is close to its end ac- 
counts for the fact that at this time, ionizations (and 
hence recombinations) take place in dense, partially neu- 
tral gas (Lyman-limit systems) and the photons pro- 
duced after recombinations are consumed inside this sys- 
* ) ionizing the IGM (see: eg. 



Furlanetto et al. 


2006 


The traction oi 


Ly-a 



2008 



1964 



usin g a combinatio n of fits tabulated by 



and 



Martin 



po et al. 
Pengelly 



(1988) for T K > 10 6 and T K < 10 3 



respectively: 

f ICC = 0.686-0.106 log 10 (T K /10 4 K)-0.009(T if /10 4 K)-°- 4 

(37) 

The luminosity density (per comoving volume) in 
Ly-a from hydrogen recombinations in the IGM, £^. M , 
is then given by 



^rec V / — Jrec' t rec- Cy Lya- 



(38) 



3.2. Lyman-a emission from excitations in the IGM 

The UV radiation that escapes the galaxies with- 
out producing ionization ends up ionizing and excit- 
ing the neutral hydrogen in the IGM and heating the 
gas around the galaxies. The high energetic elec- 
tron released after the first ionization spends its en- 
ergy in collisions/excitations, io nizations and heating the 
IGM gas until it thermalizes (Shull & van Steenberg 
1985[ ). We estimated the contribution of the direct col- 
lisions / excitations to the Ly-a photon budget and con- 
cluded that it is negligible. 

The Ly-a luminosity density due to the collisional 
emission (radiative cooling in the IGM), is given 

by: 

€c M ( z ) = n e n m q hya E Lya , (39) 

where nni = — Xi) /^ 3 ^y is the neutral hydrogen 
density, Xi is the IGM ionized fraction and QLya is the ef- 
fective collisional excitation coefficient f or Ly-a emission 
which we calculated in the same way as [Cantalup o et aL] 
( 2008 ) , but using different values for the gas temperature 
and iGM ionization fraction. 

Considering excitation processes up to the level n = 3 
that could eventually produce Ly-a emission, the effec- 
tive collisional excitation coefficient is given by: 

9Ly Q = ?l,2p + <?l,2s + Ql,3p- (40) 

The collisional excitation coefficient for the transition 
from the ground level (1) to the level {nl) is given by 

8.629 x 10~ 6 

1\,nl — 



T, 



1/2 



^Me^/^cmV 1 , (41) 



K 



combination, f rec , is temperature dependent so we used 



where f2(l, nl) is the temperature dependent effective col- 
lision strength, oj% is the statistical weight of the ground 
state, Ei jU is the energy difference between the ground 
and the nl level and ks is the Boltzmann constant. 

3.3. Scattering of Lyman-n photons emitted from 
galaxies 

Continuum emission of photons, by stars, from Ly-a to 
the Lyman-limit travels until it reaches one of the 
Ly-n lines where it gets scattered by neutral hydrogen. 
Most of this scattering will have as end result the pro- 
duction of Ly-a photons which eventually redshift out 
of the line. Since a considerable fraction of this pho- 
tons only reach a given Ly-n frequency in the IGM this 
Ly-a emission is formed as a flux that decays with r 2 
around the star that emitted the continuum photons so it 
appears diluted in frequency in line observations of point 
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sources ( Chen fc Miralda-Escude|2008 1 . This continuum 
photons are much less likely to be absorbed by the dust 
in the ISM than photons originated in recombinations. 

In intensity mapping the frequency band observed is 
much larger than in line observations so in principle all 
the continuum Ly-a photons can be detected. Using 
the Spectr al Energy Distri bution (SED) made with the 
code from Maraston ( 2005 ) we estimated that the num- 
ber of photons emitted by stars between the Ly-a plus 
the lyman alpha equivalent width and the Lyman-limit 
is equivalent to Q^m = g 31 x iq60 Mq 1 s~ x . The higher 
frequency photons are absorbed by hydrogen atoms as 
they reach the Lyman beta frequency, reemitted and 
suffer multiple scattering until they reach the Ly-a line. 
The fraction of the continuum photons emitted close to 
the Ly-a line have already redshifted to lower frequencies 
before reaching the IGM so they will not be scattered by 
the neutral hydrogen in the IGM and will not contribute 
to the radiative coupling of the 21 cm signal (they are 
already included in the calculation of the Ly-a emission 
from galaxies). 

The intensity of this emission was calculated with a 
stellar emissivity that evolves with frequency as v~ a with 
a = 0.86 and normalized to Ql™. The Ly-a luminosity 
density originated from continuum stellar radiation and 
emitted in the IGM, ^cmf, ^ s then approximately given 
by: 

« Qi G ^X yQ SFRD(z), (42) 

where the SFRD is in units of Mq per second. Note 
that in section |4j this calculation is done through a full 
simulation. 

3.4. Lyman-a Intensity 

We calculated the intensities for the several 
Ly-a sources in the IGM from their luminosity densities 
using: 

£ lGM (z) 



'IGM 



-y(z)D 



(43) 



The luminosity and so the intensity of Ly-a emission in 
the IGM depends on local values of the hydrogen ionized 
fraction, the gas temperature and the gas density. These 
parameters are correlated with each other and so theoret- 
ical calculations of the average intensity made with the 
average of this parameters may be misleading. Since this 
emission is dominated by overdense regions a clumping 
factor of a few units is usually assumed in theoretical cal- 
culations. However we decided to estimate this intensity 
without using a clumping factor since its effect can be 
extrapolated from the intensity without clumping. The 
intensity of Ly-a emission due to recombinations or col- 
lisions in the IGM is shown in figure [7] as a function of 
the hydrogen ionized fraction for different values of the 
gas temperature. 

Even for a fixed average IGM ionized fraction, the in- 
tensity of Ly-a emission is the result of emission from 
several regions and so all the values shown in figure 
are relevant. As can be seen in figure [71 the intensity 
Lyman alpha due to recombinations and collisions in the 
IGM is very sensitive to the gas temperature and to the 
fluctuations in the IGM ionized fraction. Numerical sim- 
ulations predict that the temperatures in the hydrogen 



10- 
10- 

10" 1 ' 



7 ° 10" 11 

O) 

CD 10 

3 10" 12 
> 



10" 



10- 



T K = 20000 K: exc 
T K = 20000 K: rec 
T K = 10000 K: exc 
T K = 10000 K: rec 



0.01 




1.00 



Fig. 7. — Intensity of Ly-a emission at redshift 7 due to recom- 
binations and excitations in the IGM as a function of the hydrogen 
ionized fraction. The green and red lines assume a constant gas 
temperature of 20000 K and 10000 K respectively. 



gas in the IGM can vary in the range 5000 K to 20000 K 
pave et al.|[200l] [Smith et al ||2011[ ). 

The theoretical intensities of Ly-a emission in galaxies 
and in the IGM shown in figure [^indicate that unless the 
IGM clumping factor is very high, or the Ly-a photon 
escape fraction is very low, the Lyman alpha intensity 
from the IGM at z = 7 is lower than the emission from 
galaxies. 
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Fig. 8. — Intensity of Ly-a emission at redshift seven from the 
IGM and from Galaxies as a function of the hydrogen ionized frac- 
tion and including all contributions. The green and red lines are 
the intensity of Ly-a emission in the IGM assuming a constant gas 
temperature of 20000 K and 10000 K respectively. The blue solid 
line is the intensity of Ly-a emission from Galaxies as calculated in 
the previous section. The yellow dashed lines show the intensity in 
galaxies assuming an error in A(z) of 20% due to the uncertainty 
in the ionizing photons escape fraction and due to the uncertainty 
in the emissivity of ionizing photons. The intensities in the IGM 
were calculated assuming a clu mping factor of 6 c ompatible with 
current conservative estimates. (;Pawlik et ah 2 010| . 



At higher redshifts the SFRD will decrease causing the 
Ly-a intensities from galaxies and from the IGM to de- 
crease. The escape fraction of UV photons from galaxies 
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increases as the redshift increases which will contribute 
negatively to the intensity of emission in galaxies and 
positively to the intensity of emission in the IGM. At 
high redshifts the IGM ionized fraction is small which 
contributes to a strong decrease in the intensity of emis- 
sion from the IGM compared to the intensity at z = 7. 

4. LYMAN-o INTENSITY AND POWER SPECTRUM USING 
NUMERICAL SIMULATIONS 

The intensity of Lyman alpha emission in the IGM at a 
given time and a given region is proportional to the ion- 
ized fraction, the gas temperature and the matter density 
in that region. Since these three quantities are corre- 
lated, the use of average values in the Ly-a intensity cal- 
culation highly underestimates the emission in the more 
overdense regions. Also the evolution of the average of 
the IGM ionized fraction is poorly known during the 
Epoch of Reionization (EOR). 

Some of these problems can be resolved using a com- 
putational code able to produce simulations of the IGM 
ionized fraction, the gas temperature and the matter den- 
sity in a volume high enough to properly represent our 
Universe. The use of simulations has an additional ad- 
vantage of allowing the calculation of the 3 dimensional 
power spectra of Lyman alpha emission in the IGM with- 
out the need for assuming a bias relation with the un- 
derlying dark matter distribution. 

In this section we will estimate the inhomogeneous 
Ly-a intensity from Galaxies and the IG M using a mod - 
ified version of the SimFast21 code ( Santos et al.| [2010). 
Given a set of astrophysical and cosmological parameters, 
this code is able to consistently produce 3 dimensional 
simulations of the dark matter density field, the ioniza- 
tion field, the SFRD, the scattering of Ly-n photons in 
the IGM, the X-ray heating of the IGM and even 21 cm 
spin and brightness temperature fluctuations for the sev- 
eral redshifts of the EoR. 

A proper calculation of all the heating and cooling 
mechanisms would add a high level of complexity to this 
calculation and would require a small redshift step in 
the IGM fraction calculation so we assumed a constant 
temperature in ionized regions of 10000K. Moroever, the 
results from our calculations can be easily extrapolated 
to account for a higher temperature. For example for a 
temperature of 20000K the number of recombinations in 
the IGM would decrease by a factor of 1.7 and the num- 
ber of collisions would increase more than two orders in 
magnitude. Assuming that the clumping of the IGM is 
not very high, and so Ly-a recombination emission dom- 
inates over collisional emission during most of the EOR, 
than this higher temperature would cause a small de- 
crease in the intensity of emission in the IGM and the 
Reionization period would be less extended than what 
we predict in section |4.1| We made a few modifications 
to the SimFast21 code in order to provide a consistent 
description of the ionization history and its relations to 
the Ly-a emission, which we now describe. 

4.1. IGM Ionized fraction calculation 

In the previous version of the SimFast21 code, the IGM 
ionized fraction was computed assuming that at each red- 
shift the ionization state of a region could be estimated 
from the collapsed mass in that region assuming a linear 
relation between collapsed mass and ionizing power. So 



a given spherical region o f radius R is considered ionized 
if |Furlanetto et al.|2006| ) : 

(44) 



C M coll (R) > M tot (R), 



where M co ii is the collapsed mass which corresponds to 
the total mass in halos in that region, M to t is the total 
mass in the region and £ is an ionizing efficiency param- 
eter. This efficiency parameter tries to include all the 
ionizations and recombinations produced by a halo as a 
function of its mass but has no actual physical meaning 
although its use is somewhat justified by the large un- 
certainty in the astrophysical quantities involved in the 
determination of the relation between halo mass and ion- 
izing efficiency and in the adjustment of this parameter 
in order to reproduce a reionization history compatible 
with observations. 

In order to calculate the Ly-a field however, we need to 
include the recombinations in the IGM explicitly, as well 
as directly relate the ionization process to the emitted 
stellar radiation. We therefore modified the SimFast21 
code to include these improvements. This new method 
allows a non linear relation between collapsed mass and 
ionizing power and all the parameters involved in the cal- 
culation have values based in current astrophysical con- 
straints. Also, the size of ionized regions is now set by 
the volume at which the total ionizing emissivity of the 
sources it contains equals the number of recombinations 
so that the system is in equilibrium. For each redshift 
the implementation of this method was done with the 
following steps: 

1. A halo catalog with the mass and three dimen- 
sional spatial positions was generated using the 
same method used in the original version of the 
SimFast21 code. 

2. We calculated SFRs from the halo catalogs using 
the non- linear relations in equations [24} [25j and [26} 

3. We used equation|4]to obtain the halo ionizing rate, 
./Vion, we corrected for the presence of helium using 
Ah c and multiplied it by / csc , to account for the 
photons consumed inside the galaxies. 

N^(z, M) = A Be N ion (M)f esc (z, M). (45) 

The UV ionizing rates of the halos, were then 

put in three dimensional boxes. 



4. Three dimensional boxes with the rate of recombi- 
nations in each cell, iV^ M = V ce a x n]^ 1 , were ob- 
tained from a dark matter density simulation made 
with the SimFast21 code using equation [34] with xi 
set to one and Tk = 10 4 K. 

5. Following the same procedure as in the original ver- 
sion of the SimFast21 code we applied a series of 
top-hat filters of decreasing size (this filtering pro- 
cedure was done in Fourier space) to the ionizing 
rate and the recombination rate boxes in order to 
calculate the region ionizing rate and recombina- 
tion rate. 

6. At each filtering step of radius R we found the ion- 
ized regions (HII bubbles) by checking if the region 
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ionizing rate was equal or higher than its recom- 
bination rate. With this method HII bubbles are 
always fully ionized: 



1.0 



n& m (z,r)>n; 



IGM 



(46) 



4.2. Intensity from recombinations and collisions in the 

IGM 

The SimFast21 code was built to calculate the IGM 
ionized state assuming two types of regions: one fully 
ionized (inside the HII bubbles) and other fully neutral. 
The intensity of Ly-a emission in the IGM due to re- 
combinations is a smooth function of the IGM ionized 
fraction and is dominated by emission from fully ionized 
regions (see figure[7]) so the output os the SimFast21 code 
is good enough to estimate this intensity. 

Collisions between electrons and neutral hydrogen 
atoms can also lead to Ly-a emission, however as was 
explained in section 3.2, collisional Ly-a emission only 
occurs in partly ionized regions, mainly in the the edge 
of HII bubbles, so the estimation of this emission requires 
a more detailed description of the IGM ionized state than 
the one given by the limited resolution of semi numerical 
simulations. 

Collisions are most important in regions where the 
IGM ionized fraction is locally close to 0.5 and the tem- 
peratures are high. Since high temperature regions are 
likely to be highly ionized we can deduce with the help 
of figure [7J that Ly-a emission from recombinations is 
dominant over Ly-a emission from collisions in the IGM. 




Fig. 9. — Evolution of the IGM ionized fraction as a function of 
rcdshift for the three star formation rate halo mass parameteriza- 
tions shown in section [2. 51 

depth constraint Sim2 and Sim3 have the most likely 
reionization histories and the IGM ionized fraction evo- 
lution obtained with Siml can be seen as a lower bound. 
The intensities of Lyman alpha emission from galaxies 
at redshift seven obtained with the SimFast21 code are 
similar to the more theoretical estimates summarized in 
table! 

Intensities of Lyman alpha emission in the IGM made 
with the same code are presented in table [5] 



4.3. Intensity from the scattering of Lyman-n photons 
in the IGM 

The IGM Ly-a intensity from scattering of 
Ly-n photons emitted from galaxies can also be 
calculated using data fr om the code SimF ast21. This 
code uses Equation 10 in Santos et al. (2010) to calculate 
the spherical average of the number of Ly-a photons, 
J a , hitting a gas element per unit proper area per unit 
time per unit frequency per steradian. 
The Lyman alpha intensity originated from these 
continuum photons is given by: 



IGM 



§J a E\ jya D\ 
(1 + z 



(47) 



4.4. Results 



Using the prescriptions described in the previous sec- 
tions we ran simulations Siml, Sim2 and SimS with a 
volume of 54 3 /i~ 3 Mpc 3 and 1800 cells from redshift 14 
to redshift 6. 

The obtained IGM ionization fractions, at redshift 
seven, where Xi =0.86 for simulation Siml and Xi =1.0 
for simulations Sim2 and SimS. These values are consis- 
tent with the current most likely values for t his parame- 
ter, 0.8 < x t {z = 7) < 1.0 ( |Mitra et al.|2012j ). 

The IGM ionized fraction evolution for Sim2 and for 
SimS (see figure [9]), resulted in optical depths to reion- 
ization of 0.073 and 0.082. This optical depths are 
consistent wit h the value obtained by WMAP (r = 
0.088 ± 0.015) ( |Komatsu et al |2011[ ). The optical depth 
correspondent to Siml is 0.66 which is lower than the 
current observational constraints. Based in the optical 
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TABLE 5 

Surface brightness (in observed frequency times intensity) 
of Ly-« emission from the different sources in THE IGM AT 

Z RJ 7, Z £3 8 AND z w 10. 



The intensity values found in tables [3] a nd [5] and the 
theoretical estimations plotted in figure [8j indicate that 
for the Ly-a intensity from the IGM to reach a value 
close to the emission from galaxies at z — 7 would re- 
quire a very large absorption of Ly-a photons by dust in 
galaxies. 

The resulting power spectra of Ly-a emission in galax- 
ies and in the IGM obtained with the SimFast21 code 

7 and for z 



are presented in figure 10 for z — 7 and for z = 10. 

We repeated the Ly-a power spectra calculation for 
several redshifts during the EOR and plotted the 
Ly-a power spectra as a function of redshift for several 
k in figure [TTJ 

We calculated the intensity of Ly-a emission from 
galaxies and from the IGM (intensities are shown in fig- 
ure 12 1, and found that according with our assumptions 
and as already previously seen, the Ly-a emission from 
galaxies is dominant over the Ly-a emission from the 
IGM at least during the redshift interval from z = 6 
to z = 9. 

Since the star formation halo mass relation is not very 
constrained, we can use the results obtained with Siml 
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Fig. 12. — Ly-a Intensity from galaxies (dashed lines), from the 
IGM (solid lines) as a function of redshift from our simulation 
(red thin lines) and from the theoretical calculations (blue thick 
lines). Also shown is the total Ly-a emission from the simulation 
(dotted line). All the intensities where calculated using the star 
formation halo mass relation from Siml. The theoretical intensity 
of Ly-a emission from the IGM was calculated using the average 
IGM ionization values obtained from Siml. 
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TABLE 6 

Surface brightness (in observed frequency times intensity) 

in units of [ergs _1 cm _2 sr -1 ] of LY-a emission from Galaxies 
at z k. 7. z « 8 and z « 10 for Siml, Sim2 and SimS. 



Fig. 10. — Power spectrum of Ly-a surface brightness (i/I) 
from galaxies (thin lines) and from the IGM (thick lines) at red- 
shifts 7 (top) and 10 (bottom). The shown contributions to the 
Ly-a flux are from: excitations and collisions, recombinations, con- 
tinuum emission inside the Ly-a width (from galaxies), scattering 
of Ly-n photons (from the IGM), cooling emission in galaxies and 
total emission. 




Fig. 11. — Total power spectrum of Ly-a emission during the 
EOR as a function of redshift. 
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TABLE 7 

Surface brightness (in observed frequency times intensity) 

in UNITS OF [p-™- -1 
AT2RJ7 



cm sr X ] OF LY-O EMISSION FROM THE IGM 
! AND z ss 10 FOR Siml, Sim2 AND SimS. 
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TABLE 8 

Surface brightness (in observed frequency times intensity) 
_1 ™ _2 sr~ 1 ] of total LY-a emission at z 7, 
AND z « 10 FOR Siml, Sim2 AND Sim3. 



and SimS as the lower and upper bounds to the expected 
Ly-a intensity. 

The evolution of the Ly-a intensity from galaxies, from 
the IGM and total intensity can be seen respectively in 
tables 6][7 and [8] for simulations Siml, Sim2 and Sim3. 

A map of the total Ly-a in tens ity in galaxies and in 
the IGM is presented in figure [13] for z = 7. 
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Fig. 13. — Total Ly-a Intensity from galaxies and the IGM in 
erg s — 1 cm -2 sr —1 at redshift 7. 

5. CROSS-CORRELATION BETWEEN Ly-a AND 21-CM 
OBSERVATIONS 

Observations of the 21 cm signal from the EOR will 
suffer from contamination by foregrounds and system- 
atic effects. Since both 21 cm line emission and Ly-a line 
emission trace neutral hydrogen, these two lines are ex- 
pected to be strongly correlated. The cross correlation of 
these two lines can be used as an extra method to probe 
the evolution of the IGM ionized hydrogen fraction. In 
particular the power spectra of this cross correlation will 
have a discontinuity at a scale that is related to the aver- 
age bubble size and hence the average ionization fraction 
in the Universe. 

During the EoR, the 21 cm signal from galaxies is 
much smaller than the emission from the IGM so it 
is safe to neglect both this galaxy emission and the 
shot noise emission in the cross-correlation. Since the 
Ly-a emission from galaxies is dominating over the IGM 
for most redshifts, we can just concentrate on the Ly-a - 
galaxy/2 lcm-IGM cross-correlation when analyzing the 
cross-power spectrum. The cross correlation between the 
21 cm signal and the Ly-a line in galaxies is therefore 
given by: 



Phya,2l(z,k) = IgAL^21 



Pss 



1 



1 



-Prr 



(48) 



where li\ is the average intensity of 21 cm emission, 
P Xi s(z, k) is the cross correlation power spectra between 
the ionized field and the matter density fluctuations, 
Pss(z,k) is the power spectra of matter density fluctu- 
ations and we are assuming that the Ly-a emission is a 
biased trac er o f the underlying dark matter field. 

In figure 14 we show the cross-correlation power spec- 
trum between the total Ly-a emission and the 21 cm sig- 
nal at redshifts 7, 8, 9 and 10. For simulation Siml this 
redshifts correspond to ionizing fractions of Xi = 0.86, 
Xi = 0.56, Xi = 0.35 and x; = 0.23 for redshifts 7, 8, 
9 and 10 respectively. In figure 14 the scale at which 
Phya,2i(k) goes from negative to positive is determined 
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Fig. 14. — Cross correlation power spectrum between 
Ly-a emission and the 21 cm signal for redshifts 7, 8, 9 and 10. 
Dotted lines indicate a negative correlation and solid lines indicate 
a positive correlation. 

by the average size of the ionized regions. For small scales 
the correlation is positive since fluctuations from both 
lines should be proportional to the underlying density 
fluctuations but for large scales (small k) the correlation 
is negative since the 21 cm line and the lyman alpha line 
are characteristic of neutral gas and ionized gas respec- 
tively (and there will be an extra negative contribution 
from the ionised bubbles). 

6. OBSERVATIONS 

Current observations related to Ly-a emission are re- 
stricted to narrow-band imaging of Ly-a emitters during 
reionization and the direct detection of individual emit- 
ters. This has led to order ~ 200 secure detections at 
z > 5, but systematic uncertainties remain on the frac- 
tion that are arising at low redshifts and associated with 
[OIII]/[OII] lines, among others. Due to narrow atmo- 
spheric windows, observations in near-IR bands neces- 
sary to cover the reionization epoch are also limited to 
multiple discrete bands. In any case, existing data could 
be used for a statistical study such as the power spec- 
trum to extract properties of Ly-a emitters that remain 
below the 5<r level of individual source/line identifica- 
tions. Given that detections do exist at the bright-end 
and our predictions are consistent with the Ly-a LFs de- 
rived from observational measurements, it is likely that 
a modest improvement in existing technology and pro- 
grams will lead to an experiment with sufficient sensi- 
tivity to measure the Ly-a anisotropy power spectrum 
during reionization over a broad range of redshifts. The 
main limitation, unfortunately, is that existing ground- 
based observations are very limited to small fields of view 
with narrow-bands in the redshift. 

Note that from the ground we expect a noise (vl) of 
~ 2.5 x 10 -3 erg cm~ 2 sr" 1 (assuming we can avoid the 
OH lines, otherwise, the intensity will be ~ 1.0 x 10 _1 erg 
cm~ 2 sr^ 1 ). From space, the main contamination will be 



the zodiacal light, which will have a value ~ 5 x 10 4 



erg 



cm sr . It is possible that a dedicated experiment 
from the ground can be conceived to improve our under- 
standing of reionization through detailed Ly-a mapping 
over a broad range of redshifts using specific instruments 
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Fig. 15. — Expected error on the Ly-a clustering power spectrum 
at 2 = 7 using a space based experiment. Black solid line shows the 
clustering power spectrum for Sim3 while the dashed line includes 
the shot noise. Red vertical bars shows the error. The lower blue 
solid line shows the clustering power spectrum for Siml while the 
top blue solid line shows the same for a model similar to SimS with 
the same reionization history and optical depth (from WMAP) but 
with a SFR 3 times larger and a UV escape fraction 3 times lower, 
which will generate a Ly-a luminosity function larger than what is 
usually expected at z = 7. Circles show the expected Ho power 
spectrum from z = 0.5 that will contaminate the observation and 
the crosses gives the expected Ha signal after galaxies with a Ha 
luminosity > 1.0 X 10 40 erg s" 1 are removed. 



and filters that suppress the atmospheric contamination. 
Because of this strong atmospheric contamination, sub- 
orbital and/or orbital experiments may however offer a 
better option. The predictions we have made here can 
be used as a guide in designing such instruments and 
experimen ts. 

In hgure 15 we show the expected errors at z = 7 (cen- 



tral wavelength of 0.975 /xm) for a dedicated compact 
space-borne template instrument to study Ly-a EOR 
fluctuations. We consider a 20 cm aperture and a spec- 
trometer with resolution R = A/ A A = 200. The imaging 
will be done using a 2048x2048 HgCdTe detector array 
in order to cover in one pointing a field of view of 45x45 
arcmin with a resolution of 10 arc-second pixels on the 
sky and a spectral range from 0.85 to 1.10 /im. We took 
a survey area of 20 deg 2 and a total observation time of 
2900 hours. This example shows that Ly-a EOR science 
is well within the reach of our modest template instru- 
ment. The calculated sensitivities achieved on the deep 
fields are sufficient to detect Ly-a in broad Ak/k bins 
ranging from k = 0.01 to 10 h/Mpc in both clustering 
and Poisson fluctuations. 

Ideally the spectral resolution would match the max- 
imum k available in the angular direction; however 
higher spectral resolution requires longer integration 
times needed to realize photon noise limited sensitiv- 
ity, which tends to degrade the instrument sensitivity. 
The angular resolution does not affect surface bright- 
ness sensitivity directly, but does determine the depth 
to which lower-redshift galaxies may be masked using a 
deep ancillary continuum galaxy survey. Although the 
continuum emission from galaxies can in principle be re- 
moved by looking at the signal across the frequency di- 
rection, as explained before, contamination from other 



lines at lower redshifts does poses a problem to the de- 
tection of the Ly-a signal, in particular from the Ha line. 
The most straightforward way to remove this contami- 
nation would be by masking the pixels where these low- 
z galaxies are found, either from the observation itself 
or using another, high sensitivity, continuum observa- 
tion. For this approach, the angular resolution of the 
Ly-a experiment has to be good enough in order to have 
enough pixels left after the masking. Therefore, this in- 
strument is required to have hig her angular resolution 
than spectral resolution. Figure [15] also shows the ex- 
pected contamination from the Ha line from galaxies at 
z = 0.5 (black dots). This was calculated following the 
same approach as for t he Ly-a line and using the Ha to 
SFR relation t aken from Kennicutt et al. ( 1994 ) and Ken- 
nicutt ( 1998b ). Removing low-z galaxies down to a mass 
of ~ 6.6 10 iu Mq, corresponding to a cut in Luminosity 
L > 1.0 x 10 40 erg/s, would bring this contamination 
below the Ly-a signal (black crosses) . Usi ng the Ha Lu- 



minosity function from Geach et al. (2010) normalized to 
the SFR density at z = 0.5 we get an expected angular 
density of about 25 Ha emitters per square degree per 
band, which would mean that only ~ 0.98% of the pixels 
would be masked. 

Note that the rejection of interloping low-redshift 
galaxies requires a full treatment that is beyond the 
scope of this paper. Foreground rejection may also 
be significantly enhanced by simultaneously detecting 
additional EOR spectral features beyond Ly-a , which 
are produced by interlopers with very low probability. 
Combining these Ly-a measurements with other EOR 
observations (CO, C+ and particularly HI 21 cm) 
offers additional information on EOR star- format ion, 
metallicity, and ionization history. The possibility 
of constructing an experiment in a near-IR band to 
measure the Ly-a flux in order to cor relate it with 
the 21 cm signal was also explored by W yithe et al.| 
(2007). Although, they used simple models to estimate 
the fluctuations in each of these two lines, they also 
considered several foregrounds that will contaminate the 
observations and concluded that it is possible to remove 
enough foregrounds that the intensity of radiation 
emitted from galaxies can be constrained from the cross 
correlation. 



7. SUMMARY 

In this paper we took into account the main contribu- 
tions to Ly-a emission from recombinations, collisions, 
continuum emission in galaxies and scattering of Lyman- 
n photons to calculate the intensity of Ly-a emission 
from galaxies and from the IGM during the EOR. 

We started by theoretically calculating the intensities 
using astrophysical data from several observational re- 
sults and then implemented the calculation in a simu- 
lation using a modified version of the code SimFast21 
to obtain the spatial fluctuations of Ly-a emission. The 
simulation allowed to calculate the Ly-a emission tak- 
ing into account the spatial fluctuations of the different 
astrophysical parameters, which represents an improve- 
ment over theoretical calculations that only use the av- 
erage values. 

Our simulations showed that to achieve optical depths 
compatible with the WMAP constraints the high SFRD 
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required imply that for reasonable values of UV and 
Ly-a escape fraction the intensity of Ly-a emission from 
galaxies is dominant over the emission from the IGM. 

By testing different SFR halo mass parameterizations 
we constrained the intensity of Ly-a emission from galax- 
ies to be about (1.43 - 3.57) x 1(T 8 and (4.55 - 9.73) x 
1CP 11 erg _1 cm~ 2 sr _1 at redshift 7 and 10, respectively 
which is dominant over the intensity of Ly-a emission 
from the IGM at z = 7 (about 1.6 x 1CT 5 ') but less at 
z = 10 (1.1 x 10~ 10 erg s _1 cm~ 2 sr _1 ). Since the in- 
tensity levels we found are lower than the extragalactic 
background intensity from galaxies and so are too low 
to be detected with an experiment aiming the absolute 
background intensity, we propose an intensity mapping 
experiment which will allow to measure the Ly-a power 
spectrum. 

For reasonable astrophysical conditions the process of 
hydrogen reionization was done by UV radiation origi- 
nated in galaxies with luminosities below the high red- 
shift observational threshold. In this work we showed 
the different ways by which UV emission is connected to 
Ly-a emission and so we stress how it would be useful 
to use intensity mapping of Ly-a emission to probe the 
overall intensity of UV radiation. 

Ly-a emission can also be connected to the 21 cm sig- 
nal from the Epoch of Reionization, since the continuum 



photons above the Ly-a line that redshift to this line in 
the IGM contribute to the radiative coupling of the 21 
cm signal to the gas temperature. The cross correlation 
of the Ly-a and the 21-cm lines can be used to reduce 
systematics and foregrounds encountered with 21-cm ob- 
servations. In particular the discontinuity of the cross 
correlation power spectra will provide constrains in the 
evolution of the IGM ionized fraction. 

In previous studies we have discussed the use of CO 
molecular and CII fine-structure atomic lines to comple- 
ment 21-cm data in the attempt to probe the IGM during 
reionization. Our study shows that Ly-a intensity map- 
ping is also a viable approach to probe reionization and 
is within the experimental reach over the coming decade. 
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